Cosmological Constraints from Moments of the Thermal Sunyaev-Zel'dovich Effect 
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In this paper, we explain how moments of the thermal Sunyaev-Zel'dovich (tSZ) effect can con- 
strain both cosmological parameters and the astrophysics of the intracluster medium (ICM). As the 
tSZ signal is strongly non-Gaussian, higher moments of tSZ maps contain useful information. We 
first calculate the dependence of the tSZ moments on cosmological parameters, finding that higher 
moments scale more steeply with as and are sourced by more massive galaxy clusters. Taking ad- 
vantage of the different dependence of the variance and skewness on cosmological and astrophysical 
parameters, we construct a statistic, |(T 3 )|/(T 2 ) 1 ' 4 , which cancels much of the dependence on cos- 
mology (i.e., erg) yet remains sensitive to the astrophysics of intracluster gas (in particular, to the gas 
fraction in low-mass clusters). Constraining the ICM astrophysics using this statistic could break 
the well-known degeneracy between cosmology and gas physics in tSZ measurements, allowing for 
tight constraints on cosmological parameters. Although detailed simulations will be needed to fully 
characterize the accuracy of this technique, we provide a first application to data from the Atacama 
Cosmology Telescope and the South Pole Telescope. We estimate that a Planck-like full-sky tSZ 
map could achieve a < 1% constraint on as and a la error on the sum of the neutrino masses that 
is comparable to the existing lower bound from oscillation measurements. 



I. INTRODUCTION 

The thermal Sunyaev-Zel'dovich (tSZ) effect is a spec- 
tral distortion of the cosmic microwave background 
(CMB) caused by inverse Compton scattering of CMB 
photons off hot electrons in the intracluster medium 
(ICM) of galaxy clusters [J. The tSZ effect, which is 
largest on arcminute angular scales, has traditionally 
been studied either through observations of individual 
clusters — both pointed measurements and recent 
detections in blind surveys — or indirectly through 
its contribution to the CMB power spectrum [jj, [To[ ■ 

The goal of many direct studies is to characterize the 
masses and redshifts of a well-defined sample of clusters, 
and thereby reconstruct the high-mass end of the halo 
mass function, an important cosmological quantity that 
is sensitive to a number of parameters. These parame- 
ters include ag , the rms fluctuation amplitude on scales of 
8 hr 1 Mpc. However, one can also constrain these param- 
eters through their influence on the tSZ power spectrum, 
which is exceptionally sensitive to ag in particular [Tlj . 
In recent years, many studies have obtained competitive 
constraints on ag from the power spectrum [j| [lQ] ■ This 
approach is advantageous not only because the tSZ sig- 
nal is extremely sensitive to ag, but also because it does 
not require the measurement of individual cluster masses, 
which is a difficult procedure. Unfortunately, it does re- 
quire a precise understanding of the pressure profile of 
the ICM gas for clusters over a wide range of masses and 
redshifts. Consequently, systematic errors due to theo- 
retical uncertainty in the astrophysics of the ICM have 
remained comparable to or greater than statistical errors 
in tSZ power spectrum measurements 0, [T3| • This situ- 
ation has hindered the progress of tSZ measurements in 



providing cosmological constraints. 

In this paper, we propose a method to reduce the the- 
oretical systematic uncertainty in tSZ-derived cosmolog- 
ical constraints by combining different moments of the 
tSZ effect. Although previous studies have investigated 
tSZ statistics beyond the power spectrum fl2l - [T5j . none 
have attempted to find a tSZ observable that isolates the 
dependence on either astrophysical or cosmological pa- 
rameters. In §11, we compute the variance (second mo- 
ment) and skewness (third moment) of the tSZ signal 
for a variety of ICM astrophysics scenarios while vary- 
ing erg, which allows us to probe the mass and redshift 
dependence of each statistic, as well as to derive their de- 
pendence on cosmology. We find that the skewness is not 
only more sensitive to ag than the variance, but is also 
dominated by contributions from higher-mass clusters, 
for which the ICM astrophysics is better-constrained by 
existing observations. In §111, we use these results to find 
a particular combination of the variance and skewness — 
the "rescaled skewness" — that only depends weakly on 
the underlying cosmology while remaining sensitive to 
the ICM astrophysics prescription. We test this statistic 
using a variety of models for the ICM, which should span 
the space of reasonable theoretical possibilities. In §IV, 
we apply this method to data from the Atacama Cosmol- 
ogy Telescope(ACT) QJ, Q3 and the South Pole Tele- 
scope (SPT) [3, to demonstrate its feasibility and to 
obtain a first weak constraint on the ICM astrophysics. 
In §V, we estimate the extent to which this approach can 
increase the precision of future constraints on ag and the 
sum of the neutrino masses through tSZ measurements. 
We also highlight applications to other parameters that 
affect the tSZ signal through the mass function. 

We assume a concordance ACDM cosmology through- 
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out, with parameters taking their maximum-likelihood 
WMAP5 values (WMAP+BAO+SN) [2(j unless other- 
wise specified. All masses are quoted in units of M Q //i, 
where h = ^/(lOOkms^ 1 Mpc -1 ) and H is the Hubble 
parameter today. 



II. CALCULATING TSZ MOMENTS 

A. Background 

The tSZ effect leads to a frequency-dependent change 
in the observed CMB temperature in the direction of a 
galaxy group or cluster. Neglecting rclativistic correc- 
tions (which are relevant only for the most massive sys- 
tems) , the temperature change T at angular position 
6 with respect to the center of a cluster of mass M at 
redshift z is given by [l[ 
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where the tSZ spectral function g v = xcoth(a;/2) — 4 with 
x = hv/ksTcMB, y is the Compton-y parameter, ot is 
the Thomson scattering cross-section, m e is the electron 
mass, and P e (r) is the ICM electron pressure at location 
r with respect to the cluster center. In this work we 
only consider spherically symmetric pressure profiles with 
Pe(r) = Pe(r). Also, we calculate all observables at v = 
150 GHz, where the tSZ effect is observed as a decrement 
in the CMB temperature in the direction of a cluster (i.e., 
<?„ < 0). Note that the integral in Eq. (fTJ) is taken along 
the line of sight such that r 2 = l 2 + dA(z) 2 9 2 , where dA{z) 
is the angular diameter distance to redshift z and 6 = 
\0\ is the angular distance from the cluster center in the 
plane of the sky. Given a spherically symmetric pressure 
profile, Eq. ([T]) implies that the temperature decrement 
profile is azimuthally symmetric in the plane of the sky, 
that is, T{G,M,z) = T(6,M,z). Finally, note that the 
electron pressure is related to the thermal gas pressure 
via P th = P e (5X H + 3)/2{X H + 1) = 1.932P e , where 
Xh = 0.76 is the primordial hydrogen mass fraction. 

In order to calculate moments of the tSZ effect, we 
assume that the distribution of clusters on the sky can 
be adequately described by a Poisson distribution (and 
that contributions due to clustering are negligible, which 
is valid on sub-degree scales) [13, U 
at zero lag is then given by 



The N th moment 



(T N ) = I dz 
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where dV/dz is the comoving volume per steradian at 
redshift z, dn(M, z)/dM is the comoving number density 
of halos of mass M at redshift z (the halo mass function), 
and T(6,M,z) is given by Eq. (UJ. Our fiducial integra- 
tion limits are 0.005 < z < 4 and 5 x 10 n M G /h < M < 



5 x 10 15 M Q /h. Eq. ([2]) only includes the I-halo contribu- 
tion to the N th moment at zero lag — as noted, we have 
neglected contributions due to clustering. Using the bias 
model of [H| , we find that including the 2-halo term typ- 
ically increases (T 2 ) by only 1-2%. Moreover, the 2- and 
3-halo contributions to (T 3 ) should be even less signifi- 
cant, since the higher moments are dominated by regions 
increasingly near the center of each cluster (where \T\ 
is larger). Thus, we include only the 1-halo term in all 
calculations, as given in Eq. @. 

Note that we define M to be the virial mass, that is, 
the mass enclosed within a radius r V i r [25| : 



3M 



47rA cr (z)p cr (^) 
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(3) 



where A cr (z) = 18tt 2 + S2(Q(z) - 1) - 39(ft(z) - l) 2 , 
Q(z) = n m (l + z) 3 /(Q m (l + z) 3 + OO, and p cr (z) = 
3H 2 (z) /8ttG is the critical density at redshift z. How- 
ever, some of the pressure profiles that we use below are 
specified as a function of the spherical overdensity mass 
rather than the virial mass. Thus, when necessary, we use 
the NFW density profile [26| and the concentration-mass 
relation from [27J to convert the virial mass to a spherical 
overdensity mass Ms :C /d, where M^ c (Ms t d) is the mass 
enclosed within a sphere of radius r$, c { r &,d) such that 
the enclosed density is 5 times the critical (mean matter) 
density at redshift z. To be completely explicit, this pro- 
cedure involves solving the following non-linear equation 
for r StC (or rs,d): 



r r ^<= 4 
/ 47rr' 2 /9 NFW (r', M vir ,c vir )dr' = ■znr^ c p cr (z)S (4) 
Jo 

where c V i r = r V i r /r s is the concentration parameter (r s 
is the NFW scale radius) and one replaces the critical 
density p cr {z) with the mean matter density p m (z) in 
this equation in order to obtain rs,d instead of rs, c - After 
solving Eq. (jj) to find r,5 jC , M$^ c is simply calculated via 
Ms, c = ^nrl c p cr (z)8. 

From Eqs. ([T]) and ©, it is clear that two ingredients 
are required for the tSZ moment calculation, given a spec- 
ified cosmology: (1) the halo mass function dn(M, z)/dM 
and (2) the electron pressure profile P e (r, M, z) for halos 
of mass M at redshift z. The calculation thus naturally 
divides into a cosmology-dependent component (the mass 
function) and an ICM-dependent component (the pres- 
sure profile). Because the small-scale baryonic physics 
responsible for the details of the ICM pressure profile 
is decoupled from the large-scale physics described by 
the cosmological parameters (<78;^m, •••)) ^ is conven- 
tional to determine the pressure profile from numerical 
cosmological hydrodynamics simulations (or observations 
of galaxy clusters), which are run for a fixed cosmology. 
One can then take this pressure profile and compute its 
predictions for a different cosmology by using the halo 
mass function appropriate for that cosmology. We follow 
this approach below. 
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For all the calculations in this paper, we use the M2oo,rf 
halo mass function of [28[ with the redshift-dependent 
parameters given in their Eqs. (5)— (8); we will hereafter 
refer to this as the Tinker mass function. However, the 
tSZ moments are somewhat sensitive to the particular 
choice of mass function used in the calculation. As a test, 
we perform identical calculations using the M^qq^ halo 
mass function of [Hj]. Using our fiducial cosmology and 
the Battaglia pressure profile (see below), we find that 
the M 40 o,d mass function gives a result for (T 2 ) that is 
11% higher than that found using the M 2 oo.d mass func- 
tion, while the result for (T 3 ) is 22% higher. In general, 
higher tSZ moments are more sensitive to changes in the 
mass function, as they are dominated by progressively 
higher-mass halos that live in the exponential tail of the 
mass function. While these changes are non-negligible for 
upcoming high precision measurements, they are smaller 
than those caused by changes in the choice of pressure 
profile; hence, for the remainder of this paper, we use 
the M2oo,d Tinker mass function for all calculations. 

Our calculations include four different pressure profiles 
from [2^] , |30| . and [Til ] (and additional results using the 
profile from [34| ) , which we briefly describe in the follow- 
ing. We consider two profiles derived from the simula- 
tions of [3l|, which are reported in [29) . The first, which 
we hereafter refer to as the Battaglia profile, is derived 
from hydrodynamical simulations that include radiative 
cooling, star formation, supernova feedback, and feed- 
back from active galactic nuclei (AGN). These feedback 
processes tend to lower the gas fraction in low-mass clus- 
ters, as gas is blown out by the injection of energy into 
the ICM. The smoothed particle hydrodynamics used in 
these simulations also captures the effects of non-thermal 
pressure support due to bulk motions and turbulence. 
The second profile that we use from [2^| , which we here- 
after refer to as the Battaglia Adiabatic profiltj]], is de- 
rived from hydrodynamical simulations with all forms of 
feedback turned off. This setup leads to a non-radiative 
adiabatic cluster model with only formation shock heat- 
ing present. Thus, in these simulations, the gas fraction 
in low-mass clusters is close to the cosmological value, 
f4/f2 m (as is the case in high-mass clusters). This leads 
to much higher predicted amplitudes for the tSZ mo- 
ments (see below). Note that the Battaglia profile is 
specified as a function of M200.0 while the Battaglia Adi- 
abatic profile is given as a function of Msoo.c- 

In addition, we use the "universal" pressure profile de- 
rived in [30| . which we hereafter refer to as the Arnaud 
profile. This profile is obtained from a combination of 
XMM-Newton observations of massive, z < 0.2 clusters 
and hydrodynamical simulations that include radiative 
cooling and some feedback processes (though not AGN 
feedback). The simulations are used to extend the profile 
beyond -R500, due to the lack of X-ray photons at large 
radii. Importantly, the normalization of this profile is 



1 N. Battaglia, priv. comm. (to be included in |29|Q 



Profile 


A 2 [/iK 2 ] 


A z htK 3 ] 


a.2 


a?3 


Arnaud 


20.5 


-1790 


7.9 


11.5 


Battaglia 


22.6 


-1660 


7.7 


11.2 


Batt. Adiabatic 


47.1 


-3120 


6.6 


9.7 


Komatsu-Seljak 


53.0 


-3040 


7.5 


10.6 


Shaw 


23.8 


-1610 


7.9 


10.7 



TABLE I. Amplitudes and power-law scalings with as for the 
tSZ variance and skewness, as denned in Eq. ([5}. The first 
column lists the pressure profile used in the calculation (note 
that all calculations use the Tinker mass function). The am- 
plitudes are specified at erg = 0.817, the WMAP5 maximum- 
likelihood value. All results are computed at v = 150 GHz. 

obtained using hydrostatic equilibrium (HSE)-based es- 
timates of the mass Mmo, c , which are known to be biased 
low by w 10 — 15% [12l ■ Thus, following [HI, we use a 
HSE-bias correction of 13% in our calculations with this 
profile, that is, we set M^ff = 0.87M 50 o, c . This cor- 
rection slightly lowers the amplitude of the tSZ moments 
for this profile. 

We also use the profile derived in [ll[ , which we here- 
after refer to as the Komatsu-Seljak (or K-S) profile. This 
profile is derived analytically under some simplifying as- 
sumptions, including HSE, gas tracing dark matter in the 
outer regions of clusters, and a constant polytropic equa- 
tion of state for the ICM gas. In particular, it includes no 
feedback prescriptions or sources of non-thermal pressure 
support. The gas fraction in low-mass clusters (indeed, in 
all clusters) is thus equal to the cosmological value, lead- 
ing to high predicted amplitudes for the tSZ moments. 
This profile is now known to over-predict the tSZ power 
spectrum amplitude [H, [l(| , but we include it here as an 
extreme example of the possible ICM physics scenarios. 
It is specified as a function of the virial mass M as defined 
in Eq. ©. 

Finally, we present results computed using the fiducial 
profile from [34j , which have been graciously provided by 
the authors of that study. This profile is based on that 
of [llj], and thus we hereafter refer to this as the Shaw 
profile. It is derived from an analytic ICM model that ac- 
counts for star formation, feedback from supcrnovae and 
AGN, and non-thermal pressure support. Its tSZ predic- 
tions are generally similar to those of the Battaglia pro- 
file. In particular, it also predicts a suppression of the gas 
fraction in low-mass clusters due to feedback processes. 

Overall, then, we have two profiles that include a va- 
riety of detailed feedback prescriptions (Battaglia and 
Shaw), one profile based primarily on an empirical fit to 
X-ray data (Arnaud) , and two profiles for which the ICM 
is essentially in HSE (Battaglia Adiabatic and Komatsu- 
Seljak). The Arnaud, Battaglia, and Shaw models all 
agree reasonably well with X-ray observations of massive, 
low-redshift clusters [35|. The Komatsu-Seljak profile is 
somewhat discrepant for these high-mass systems, but 
disagrees more significantly with the predictions of the 
feedback-based profiles for low-mass groups and clusters. 
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FIG. 1. The tSZ variance versus ag obtained from five dif- 
ferent pressure profiles (using the Tinker mass function) and 
one direct simulation measurement. The scalings with erg are 
similar for all the models: (T 2 ) oc erg' 6-7,9 . It is evident that 
(T 2 ) is very sensitive to erg, but the scatter due to uncer- 
tainties in the ICM astrophysics (which is greater than twice 
the signal for large og) makes precise constraints from this 
quantity difficult. 



The predicted gas fraction in these low-mass objects is 
a major source (along with non-thermal pressure sup- 
port) of the difference in the tSZ predictions between the 
feedback- or X-ray-based models (Arnaud, Battaglia, and 
Shaw) and the adiabatic models (Battaglia Adiabatic and 
Komatsu-Sclj ak) . 

We also analyze a tSZ simulation from [36| that cov- 
ers an octant of the sky, providing a nontrivial test of 
our calculations. We hereafter refer to this as the Se- 
hgal simulation. The simulated data was produced by 
populating a large dark matter iV-body simulation with 
gas according to a polytropic cluster model that also in- 
cludes some feedback prescriptions. However, the model 
requires all clusters — including low-mass objects — to 
contain enough gas to reach the cosmological value of 
the gas fraction. Note that the simulation results include 
tSZ signal from components that are not accounted for 
in the halo model-based calculations, including substruc- 
ture within halos, deviations from the globally averaged 
pressure profile, and diffuse emission from the intergalac- 
tic medium. These effects are expected to contribute to 
the tSZ power spectrum at the w 10 — 20% level at high- 
l . The simulation was run using <j% = 0.8, with the 
other ACDM parameters taking values consistent with 
WMAP5. Since an average pressure profile has not been 



FIG. 2. Similar to Fig. 1, but now showing the tSZ skewness. 
The scalings with erg are again similar for all the models: 
(T 3 ) oc erg' 7-11,5 . As for the variance, the sensitivity to erg is 
pronounced, but degraded due to uncertainties in the ICM as- 
trophysics (as represented by the different choices of pressure 
profile) . 



derived from this simulation, we cannot re-compute its 
predictions for different cosmologies, and thus it is pre- 
sented as a single data point in the figures throughout 
this paper. 

B. Results 

We compute Eq. ([2]) for the variance (N = 2) and the 
unnormalized skewness (N = 3) for each of the profiles 
while varying as and keeping the other cosmological pa- 
rameters fixed. The results are shown in Figs. 1 and 
2. (See Appendix B for results involving the unnormal- 
ized kurtosis, i.e., TV = 4.) We find that the scalings of 
the variance and skewness with as are well-described by 
power-laws for each of these profiles, 




where we have normalized to the WMAP5 value of cr 8 . 
The amplitudes ^2,3 and scalings 012.3 for each of the 
pressure profiles are given in Table I. The scalings are 
similar for all the pressure profiles. Note that the slightly 
steeper scalings for the Battaglia and Shaw profiles are 
due to their inclusion of feedback processes that suppress 
the tSZ signal from low-mass clusters (the Arnaud pro- 
file also suppresses the signal from these objects, though 
it is primarily based on a fit to higher-mass clusters). 



Battaglia - - K-S M -«* [M ° /h] Battaglia - - K-S M -«* [M s /h] 

Arnaud - ■ Shaw Batt. Adia. Arnaud - ■ Shaw Batt. Adia. 



FIG. 3. Fraction of the tSZ variance contributed by clusters 
with virial mass M < M max . 



Thus, (T 2 ) and (T 3 ) for these profiles are dominated by 
higher-mass (rarer) objects than for the Komatsu-Scljak 
or Battaglia Adiabatic profiles, and they are correspond- 
ingly more sensitive to erg. Despite the steep scalings, it 
is clear in Figs. 1 and 2 that the systematic uncertainty 
due to the unknown ICM astrophysics significantly de- 
grades any potential constraints that could be derived 
from these observables. We address possible ways around 
this problem in the next section. 

We also investigate the characteristic mass and red- 
shift ranges responsible for the variance and skewncss 
signals. Figs. 3 and 4 show the fraction of the vari- 
ance and skewness signals contributed by clusters with 
virial mass M < M max . We find that the variance re- 
ceives ~ 50 — 60% of its amplitude from clusters with 
M < 2-3 x 10 14 M Q //i, while the skewness receives only 
« 20 — 40% of its amplitude from these less massive ob- 
jects. These results vary for the different profiles, as the 
different feedback prescriptions suppress the amplitude 
contributed by low-mass clusters by different amounts. 
More massive clusters are dominated by gravitational 
heating and are thus less sensitive to energy input via 
feedback from AGN, turbulence, and other mechanisms 
d^H^l- Since the Komatsu-Seljak and Battaglia Adia- 
batic profiles include no feedback, the amplitudes of their 
variance and skewness are somewhat more weighted to- 
ward low-mass objects. In all cases, though, these results 
indicate that the ICM astrophysics underlying the tSZ 
skewness is much better constrained by observation than 
that responsible for much of the variance (or power spec- 
trum), as the skewness amplitude is dominated by sig- 



FIG. 4. Fraction of the tSZ skewness contributed by clusters 
with virial mass M < iW max . Comparison with Fig. 3 indi- 
cates that the skewness signal arises from higher-mass clusters 
than those that comprise the variance. 



nificantly more massive objects which have been better 
studied observationally. This explains the smaller scatter 
seen between the different profiles in Fig. 2 as compared 
to Fig. 1: the skewness signal is dominated by objects for 
which the ICM astrophysics is less uncertain. 

In addition, we find that the variance receives « 25% 
(Arnaud/Battaglia) - 45% (Battaglia Adiabatic) of its 
amplitude from groups and clusters at z > 1, while 
the skewness receives only ss 15% (Arnaud/Battaglia) 
- 35% (Battaglia Adiabatic) of its amplitude from these 
high-redshift objects. (The Komatsu-Seljak results lie in 
the middle of these ranges.) These results imply that a 
greater fraction of the tSZ skewness than the variance 
originates from clusters that have been studied using X- 
rays, lensing, and other techniques. On the other hand, 
many of the objects comprising the variance (or power 
spectrum) have not been directly observed. 

As it depends sensitively on ag, and is less affected by 
uncertainties in astrophysical modeling than the power 
spectrum, the tSZ skewness is a powerful cosmologi- 
cal probe. Nevertheless, despite the signal originating 
from characteristically higher-mass, lower- redshift clus- 
ters, the amplitude of the skewness is still quite uncertain 
due to the poorly understood astrophysics of the ICM, 
as can be seen in Fig. 2. We investigate methods to solve 
this problem in the remainder of this paper. 
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FIG. 5. The rescaled skewness §3 t p for j3 — 1.4 plotted against 
as. It is evident that this statistic is nearly independent of 
erg, as expected based on the scalings in Table I. However, it 
is still dependent on the ICM astrophysics, as represented by 
the different pressure profiles. Thus, this statistic can be used 
for determining the correct gas physics model. 



III. ISOLATING THE DEPENDENCE ON ICM 
ASTROPHYSICS 



FIG. 6. Fraction of the rescaled skewness §3 t p for /3 = 1.4 
contributed by clusters with virial mass M < M max . A sig- 
nificant fraction of the signal comes from low-mass objects 
— even more so than the variance (Fig. 3). One can thus 
interpret this statistic as a measure of the gas fraction in low- 
mass clusters: pressure profiles that yield f gas ~ fij,/fl m in 
low-mass objects lead to a small value for this statistic (since 
they give much larger values of (T )), while pressure profiles 
that include significant feedback effects — thus lowering the 
gas fraction in low-mass objects — lead to larger values for 
this statistic. 



As we have described, the cosmological utility of tSZ 
measurements is limited by their sensitivity to unknown, 
non-linear astrophysical processes in the ICM, despite 
their high sensitivity to the underlying cosmology. How- 
ever, this theoretical uncertainty can be minimized by 
combining measurements of the tSZ skewness with mea- 
surements of the tSZ power spectrum or variance. The 
use of multiple probes provides an additional statistical 
handle on both ICM astrophysics and cosmology. 

If the tSZ variance and skewness depend differently 
on the ICM astrophysics and background cosmology, it 
may be possible to construct a statistic that "cancels" 
the dependence on one or the other. Here, we focus on a 
statistic that cancels the dependence on cosmology, but 
preserves a dependence on the ICM astrophysics. For 
an alternative approach that attempts to directly cancel 
the ICM astrophysics dependence, see Appendix A. In 
particular, we consider the following statistic, which we 
call the "rescaled skewness" : 



S: 



3,/3 



KT3)| 

(T 2) 



2\/3 



(0) 



where the value of f3 is left to be determined. From the 
results in Table I, we know that (T 2 ) oc <jg' 6 ~ 7 ' 9 and 



(T 3 } 



oc a 



9.7-11.5 



for these pressure profiles. Thus, if we 



choose ft = 1.4, then the resulting statistic will cancel 
the dependence on a&, and thus should be nearly inde- 
pendent of the background cosmology, since as is the 
dominant cosmological parameter for the tSZ observ- 
ables. Nonetheless, some dependence on the ICM as- 
trophysics may remain. We investigate this dependence 
in Fig. 5. It is evident that the rescaled skewness with 
p = 1.4 is nearly independent of as, as expected based on 
the scalings. However, this statistic still shows a depen- 
dence on the gas physics model — for eg = 0.817, the 
scatter in Ss^ = ia between the different pressure pro- 
files is « 35%. For the same value of erg, the scatter 
in (T 2 ) is » 50% between these profiles, while for (T 3 ) 
it is « 35%. Thus, the scatter in this statistic is es- 
sentially identical to that in the skewness. The rescaled 
skewness with j3 = 1.4 hence has a very useful property: 
because it only depends weakly on erg, it can be used to 
determine the correct ICM gas physics model, nearly in- 
dependent of the background cosmology. Although the 
value of (3 could be chosen slightly differently to mini- 
mize the small residual dependence on eg for any of the 
individual pressure profiles, we find that (3 = 1.4 is the 
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best modcl-indcpendcnt choice available, especially given 
current levels of observational precision and theoretical 
uncertainty. After fixing the value of /?, one can proceed 
to measure this statistic from the data (see §IV). 

Moreover, the ordering of the results in Fig. 5 suggests 
a possible interpretation of this statistic. The pressure 
profiles for which f gas w f2(,/O m in all halos — includ- 
ing low-mass groups and clusters — yield low values for 
5 , 3,/9=i.4 J while the profiles with feedback prescriptions 
that lower f gas in low-mass objects yield high values of 
this statistic. This can be explained by the fact that the 
additional contributions from low-mass clusters lead to 
higher values of (T 2 ) in the former set of profiles, while 
(T 3 ) is not significantly affected, since it is dominated 
by contributions from more massive objects. Looking 
at Eq. ©, it then follows that the profiles without f gas 
suppression have lower values for £3, ,8=1.4. Thus, the 
rescalcd skewness with /3 = 1.4 is a measure of the typi- 
cal gas fraction in low-mass groups and clusters. 

We verify this interpretation by calculating the char- 
acteristic mass range responsible for the S3 ,9=1.4 signal. 
Fig. 6 shows the fraction of this signal contributed by 
clusters with virial mass M < M max . As anticipated, it 
receives significant contributions from very low-mass ob- 
jects: w 30 — 50% of the amplitude comes from clusters 
with M < 6-7 x 10 13 M Q /h (this depends quite strongly 
on the particular profile used, as seen in the figure). A 
measurement of this quantity will thus constrain the av- 
erage amount of gas in these low-mass, high-redshift ob- 
jects, which have not yet been observed by other tech- 
niques. We present a first application of this technique 
to data from ACT and SPT in the following section. 



IV. APPLICATION TO ACT AND SPT DATA 

ACT and SPT have measured the tSZ signal in maps 
of the microwave sky with arcminute angular resolution. 
In order to apply the rescaled skewness with (3 = 1.4 to 
the data, we rephrase this statistic to match the observa- 
tional quantities as closely as possible. In particular, we 
use the SPT measurement of C3000 03 , the amplitude of 
the tSZ power spectrum at £ = 3000, and the ACT mea- 
surement of (T 3 ) [13, the filtered tSZ skewness (defined 
more precisely below). 

We use C3000 rather than (T 2 ) in order to circumvent 
uncertainties that arise when converting between the tSZ 
power spectrum and variance. ACT and SPT effectively 
measure the amplitude of the tSZ power spectrum only 
at scales of order £ = 3000, and therefore a template for 
Ci must be used in order to calculate the tSZ variance 
from this measurement. Such a template relies on an 
ICM physics model, and since this model is what we are 
hoping to constrain, we circumvent the calculation of the 
variance by working directly with C3000 instead of (T 2 ). 
We use the SPT measurement of C3000 due to its higher 
signal-to-noise ratio (SNR) than the measurement from 
ACT. Finally, note that SPT measures this quantity at 
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FIG. 7. Similar to Fig. 5, but calculated in terms of observed 
tSZ statistics, namely, C3000 (the amplitude of the tSZ power 
spectrum at £ — 3000) and (T 3 ) (the filtered skewness as de- 
fined in [HI]). The green point shows the current constraint 
using the SPT measurement of C3000 QU and the ACT mea- 
surement of (f 3 ) [3i|. 



an effective frequency of 152.9 GHz [T(|, and thus we 
perform the relevant calculations at this frequency. 

ACT recently reported the first detection of a higher- 
point tSZ observable: the filtered skewness, (T 3 ) 
This quantity is similar to the tSZ skewness that we 
calculate above, but its value has been modified by an 
^-space filter applied to the ACT maps, as well as a tem- 
perature fluctuation cut-off used to remove outlying pix- 
els in the ACT maps. We explicitly account for these 
steps in our calculations. First, we Fourier transform the 
projected temperature decrement profile in Eq. (UJ and 
apply the £-space filter used in [38[ to each "cluster" of 
mass M and redshift z in the integrals of Eq. (|2|). We 
then inverse Fourier transform to obtain the filtered tem- 
perature decrement profile for each cluster in real space. 
Second, we place each cluster in an idealized ACT pixel 
and compute the observed temperature decrement, ac- 
counting carefully for geometric effects that can arise de- 
pending on the alignment of the cluster and pixel centers. 
If the observed temperature decrement exceeds the 12a 
cut-off used in [38[ , then we discard this cluster from the 
integrals. We thus replicate the data analysis procedure 
used in the ACT measurement as closely as possible. The 
net effect is to reduce the value of the tSZ skewness by 
up to 90 — 95%. The reduction comes almost entirely 
from the £-space filtering; the cut-off used in the second 
step only has a small effect. Finally, note that the ACT 
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FIG. 8. The tSZ variance versus the sum of the neutrino 
masses Em„, with A^ = 2.46 x 1(T 9 (its WMAP5 value). 
The dependence is not precisely captured by a simple scaling 
as for ag (note that the axes are log- linear), but the curves 
are well-fit by quadratic polynomials. 



FIG. 9. The tSZ skewness versus the sum of the neutrino 
masses Em„, with A^ = 2.46 x 10" 9 (its WMAP5 value). 
The dependence is not precisely captured by a simple scaling 
as for erg (note that the axes are log-linear), but the curves 
are well- fit by cubic polynomials. 



measurement is at an effective frequency of 148 GHz [9|, 
and thus we perform the relevant calculations at this fre- 
quency. 

We use these calculations to construct a statistic anal- 
ogous to S , 3 i( 9=i.4, but defined in terms of the observa- 
tional quantities: \(T 3 )\/C^ . In order for this statistic 
to work in the manner seen in Fig. 5, we thus require the 
scalings of C3000 and (T 3 ) with cr 8 to be close to those 
reported in Table I for (T 2 ) and (T 3 ), respectively. Our 
calculations verify this claim, as can be seen immediately 
in Fig. 7. 

Fig. 7 shows the result from ACT and SPT plotted 
against the results of our theoretical calculations. Un- 
fortunately, the error bar on the data point is too large 
to deduce a preference for any particular pressure profile 
(though the central value is closer to the profiles with sig- 
nificant feedback than those with an unsuppressed value 
of fgas in all halos). Note that the error bar includes sig- 
nificant contributions due to cosmic variance (i.e., sam- 
ple variance) resulting from the limited sky coverage of 
ACT and SPT, as detailed in [H and [To|. The nearly 
full-sky results from Planck will have far smaller cosmic 
variance contributions to the error. Moreover, the error 
is currently dominated by the uncertainty on (T 3 ), which 
will greatly decrease in the near future — an upcoming 
SPT measurement should increase the SNR on (T 3 ) by 
a factor of 3 (3~4| . 

Nonetheless, other theoretical issues must still be over- 



come in order to fully characterize this technique. In par- 
ticular, the tSZ power spectrum receives non-negligible 
contributions (up to 15% at £ = 3000) due to deviations 
about the mean pressure profile found in simulations [29| . 
The tSZ skewness likely receives similar contributions, 
which are not accounted for in our halo model-based cal- 
culations. Future studies incorporating the results of hy- 
drodynamical cosmological simulations will be necessary 
to fully understand this method. Fig. 7 is simply a proof 
of concept for this technique, which may soon give inter- 
esting constraints on the ICM astrophysics of low-mass 
groups and clusters. By using the constrained pressure 
profile to interpret measurements of the power spectrum 
or skewness, it will be possible to derive stronger con- 
straints on eg from tSZ measurements. 

V. FUTURE COSMOLOGICAL CONSTRAINTS 

A. Extension to other parameters 

While we have focused on erg, the techniques that we 
have outlined are in principle sensitive to many cosmo- 
logical parameters. These include the standard parame- 
ters £lbh 2 , to which the tSZ signal is somewhat sensitive 
(though this parameter can be very well constrained by 
the primordial CMB), and f2 m , to which the tSZ sig- 
nal is not particularly sensitive (especially in the range 
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FIG. 10. The tSZ variance versus f NL , with A\ = 2.46 xlCT 9 
(its WMAP5 value). The dependence is not precisely cap- 
tured by a simple scaling as for erg (note that the axes are 
log- linear), but the curves are well- fit by quadratic polynomi- 
als. 
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FIG. 11. The tSZ skewness versus f NL , with A^ = 2.46 x 
10 -9 (its WMAP5 value). The dependence is not precisely 
captured by a simple scaling as for <Tg (note that the axes are 
log- linear), but the curves are well-fit by cubic polynomials. 



0.15 < Q m < 0.4) [ll|. Here, we instead focus on cur- 
rently unknown parameters, to which the tSZ moments 
are sensitive through the halo mass function. Such pa- 
rameters include the sum of the neutrino masses Em„, 
the non-Gaussianity parameters (e.g., Jnl), and the dark 
energy equation of state w(z). All of these parameters 
modify the number of massive clusters in the low-redshift 
universe: for example, non-zero neutrino masses suppress 
the number of clusters, while positive Jnl provides an 
enhancement. Thus, these parameters also suppress or 
enhance the amplitude of the tSZ moments. 

As a first example, we compute the tSZ variance 
and skewness for a fixed WMAP5 background cosmol- 
ogy with massive neutrinos added. To be clear, we fix 
= 2.46 x 10~ 9 (its WMAP5 value), not a s = 0.817, 
as the presence of massive neutrinos will lead to a lower 
value of as inferred at the low redshifts from which the 
tSZ signal originates. We compute the change in the 
mass function due to the neutrinos following a prescrip- 
tion similar to that of [39j , except that we input the sup- 
pressed linear theory matter power spectrum to the Tin- 
ker mass function rather than that of as was done 
in [3^] . A more precise recipe can be found in j4j [42j , 
but this approach should capture the relevant physical 
effects. 

The results of this analysis are shown in Figs. 8 and 9. 
As expected, the tSZ variance and skewness both exhibit 
sensitivity to Hm v , although it is currently overwhelmed 
by the uncertainty in the ICM astrophysics. However, 



using the method outlined in the previous section to con- 
strain the ICM astrophysics and thus break its degener- 
acy with cosmology, a measurement of the sum of the 
neutrino masses with the tSZ effect might be feasible. 

As a second example, we compute the tSZ variance and 
skewness for a fixed WMAP5 background cosmology with 
a non-zero value of /jvl added. Again, to be clear, we fix 
A|, = 2.46 x 10~ 9 (its WMAP5 value), not ct 8 = 0.817, 
as the non-zero value of /jvl will change the value of erg 
inferred at the low redshifts from which the tSZ signal 
originates. We compute the change in the mass function 
due to primordial non-Gaussianity using the results of 
[43j. Although this is technically a "friends-of- friends" 
(FOF) mass function, we follow [44[ in assuming that the 
ratio of the non-Gaussian mass function to the Gaussian 
mass function is nearly universal (and hence applicable 
to both FOF and spherical overdensity mass functions), 
even if the underlying mass function itself is not. While 
this assumption may not be valid at the percent level, 
our approach should capture the relevant physical effects 
of primordial non-Gaussianity on the tSZ moments. In 
particular, Jnl > leads to a significant increase in the 
number of massive halos at late times, while /ml < has 
the opposite effect [43|, H||. Thus, the amplitudes of the 
tSZ statistics should increase or decrease accordingly, as 
they are sourced by these massive halos. 

The results of this analysis are shown in Figs. 10 and 
11. As expected, the tSZ variance and skewness both 
exhibit sensitivity to Jnl, although for the currently al- 



10 



lowed range of values (—10 < Jnl < 74 [46| for the "lo- 
cal" shape) it is overwhelmed by the uncertainty in the 
ICM astrophysics. However, using the method outlined 
in the previous section to constrain the ICM astrophysics 
and thus break its degeneracy with cosmology, a tight 
constraint on primordial non-Gaussianity with the tSZ 
effect might be achievable. In addition, it is worth noting 
that the tSZ moments are sensitive to non-Gaussianity 
on cluster scales, rather than the large scales probed 
by the primordial CMB or large-scale halo bias. Thus, 
the tSZ signal can constrain /jyx on much smaller scales 
than other observables, allowing for complementary con- 
straints that are relevant for models of inflation that pre- 
dict strongly scale-dependent non-Gaussianity [47}. 

B. Estimating the constraints from Planck 

We briefly estimate the expected constraints on erg and 
YiTTiy using tSZ measurements from Planck. Planck's 
spectral coverage will allow the separation of the tSZ sig- 
nal from other CMB components, likely yielding a full- 
sky tSZ map. However, given issues with bandpass un- 
certainties and CO contamination, it is difficult to predict 
the SNR for a Planck measurement of S3 p=xA- Based 
on the Sehgal simulation analysis above (which covers an 
octant of the sky) , we anticipate that a cosmic variance- 
limited full-sky measurement could achieve a SNR 90 
for the tSZ variance or a SNR « 35 for the tSZ skewness. 
For 5 l 3 i/ 3 = i.4, we estimate a SNR ea 55, which implies that 
the ICM astrophysics can be constrained to the < 5% 
level. This result will allow for a < 1% constraint on erg 
after subsequently applying the correct model to inter- 
pret tSZ measurements. Note that these estimates fully 
account for both Poisson and cosmic variance error, as 
they are derived from a cosmological simulation. Techni- 
cally, we have neglected the fact that the cosmic variance 
error will scale with eg (see the discussion in [Hj]), but 
as long as the true value is not so far from erg = 0.8 as to 
conflict with all recent measurements of this parameter, 
our estimates will be accurate. 

Forecasting a £m„ constraint is somewhat more diffi- 
cult because it is highly degenerate with <jg, and we have 
not accounted for other possible parameter dependences 
(e.g., one might expect fl m h 2 to enter through its role 
in the mass function). However, as described in [39|, the 
effect of massive neutrinos is essentially captured by a 
low-redshift suppression of erg as compared to the value 
inferred from the primordial CMB at z = 1100 (i.e., A^,). 
This line of reasoning was also used in the neutrino mass 
constraint derived by the Chandra Cluster Cosmology 
Project [48[, who found an upper bound of Em^ < 0.33 
cV at 95% CL. Based on the results of [H, Sm„ « 0.1 
cV is roughly equivalent to a w 3% decrease in er 8 as 
inferred at z = 1 compared to z = 1100. The WMAP5 
measurement of erg has an uncertainty of as 3% [2(| , and 
the equivalent result from Planck should be < 1% [49| . 
If the technique that we have outlined works as planned, 



it could yield a < 1% measurement of erg a t low red- 
shift. This corresponds to a sensitivity to Y>m v ps 0.1-0.2 
cV. Given the known lower bound of £m„ > 0.05 eV, it 
might even be possible to detect the sum of the neutrino 
masses through tSZ measurements. Clearly the Planck 
results will have a smaller SNR than the full-sky, cosmic 
variance- limited assumption made above — and there are 
important degeneracies between the relevant parameters 
— but nonetheless these methods show great promise in 
constraining any parameter that affects the tSZ signal 
through the mass function. 
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Appendix A: Directly Canceling the Gas Physics 

In this appendix, we present a different approach to 
circumventing the degeneracy between ICM astrophysics 
and cosmological parameters in tSZ measurements. In 
particular, we investigate the rescaled skewness statistic 
in Eq. © for other values of the exponent (3 to see if 
they might directly cancel the gas physics dependence, 
but preserve a sensitivity to the underlying cosmology. 
We heuristically motivate the magnitude of this exponent 
in the following way. We consider the effect of increasing 
the amount of energy input through feedback on both the 
tSZ variance and skewness in a theoretical model for the 
ICM. Increasing the feedback should substantially reduce 
the signal from low-mass clusters, and generally affect the 
pressure profile in the outer regions of clusters. However, 
high-mass clusters are much less affected by feedback, so 
their tSZ signal should be essentially unchanged. This 
fact implies that the variance, which depends more on 
the signal from low-mass clusters (as shown earlier) and 
receives contributions at large distances from the clus- 
ter center, should be reduced much more by feedback 
than the skewness, which depends more on the signal 
from massive clusters and is mostly produced in the in- 
ner regions of clusters. In order to construct a rescaled 
skewness in Eq. © that is unchanged as the amount of 
astrophysical feedback is increased, we expect to expo- 
nentiate the variance (which is more sensitive to feed- 
back) to a lower power than the skewness, in order to 
compensate. Thus, we anticipate /3 < 1 for a statistic 
that is somewhat insensitive to the ICM astrophysics, 
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FIG. 12. The rescaled skewness S3 t p for = 0.7 plotted 
against as- Clearly there is significantly less scatter for this 
statistic between the different pressure profiles, but it still 
scales as 0"g -6 . Thus, this value of /3 provides a method for 
constraining cosmological parameters that is fairly indepen- 
dent of the details of the gas physics. The simulation point is 
somewhat low, but the disagreement is not statistically sig- 
nificant. This point was also not used in determining the 
optimal value of p. 



but remains sensitive to the underlying cosmology. 

To derive this exponent quantitatively, we find the 
value of (3 that minimizes the scatter between different 
choices of pressure profile for the same values of the cos- 
mological parameters. We focus specifically on erg, as 
the tSZ moments are most sensitive to this parameter. 
For 0.73 < erg < 0.9 (and all other parameters taking 
their WMAP5 values), we find that (3 = 0.7 minimizes 
the scatter between the different pressure profiles. As 
these profiles span a wide range of ICM prescriptions and 
tSZ power spectra predictions, we expect this result to 
be fairly robust. We verify this claim by computing the 
rescaled skewness with f3 = 0.7 for the Sehgal simulation. 
The results are shown in Fig. 12. Although the simula- 
tion point is somewhat low, it is only la away from the 
halo model-based results at erg = 0.8. 

In order to gain more understanding of the rescaled 
skewness with f3 — 0.7, we also investigate the charac- 
teristic mass scales that contribute to this statistic. The 
results are shown in Fig. 13. For a WMAP5 cosmology, 
we find that §3^=0.7 receives ~ 25 — 50% of its amplitude 
from clusters with M < 2-3 x 10 14 M Q /h. In general, the 
signal for S3 ,,3=0.7 is dominated by mass scales similar 
to those that comprise the skewness, i.e., higher masses 
than those responsible for the variance. It is interest- 
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FIG. 13. Contributions to the rescaled skewness with /3 — 
0.7 from clusters of various masses, for each of the different 
pressure profiles. The contributions to this statistic are very 
similar to those found for the skewness in Fig. 4, as one would 
expect based on the definition of S3, ,3=0. 7- 



ing to note that this statistic receives contributions from 
rather different mass scales for the different profiles in 
Fig. 13, even though the final result is fairly independent 
of the profile choice (as seen in Fig. 12). 

We also note that the same statistic with j3 = 0.7 be- 
haves similarly when applied to £m„, as seen in Fig. 14. 
To be clear, we do not re-fit for the value of (3 that mini- 
mizes the difference between the various pressure profiles 
in this plot. If we did perform such a fit, the best-fitting 
exponent to minimize the dispersion is again f3 = 0.7, in 
exact agreement with the best-fitting value for erg. 

Despite the promising nature of the S3. ^=0.7 statistic, 
we must emphasize that we lack a strong theoretical mo- 
tivation for this quantity, and thus one must be cautious 
in assuming its insensitivity to changes in the gas physics. 
For example, if one alters the amount of non-thermal 
pressure support by changing an overall parameter that 
affects the normalization of all halos' pressure profiles 
(e.g., the «o parameter in the Shaw model [Hj]), then 
this statistic will be affected in a nontrivial manner. It 
is likely that requiring the model to match the observed 
X-ray data for massive, low-rcdshift clusters will prevent 
a drastic change along these lines, but it is still a possibil- 
ity. In addition, as mentioned earlier, the tSZ skewness 
is fairly sensitive to the choice of mass function used in 
these calculations (though the variance is less affected), 
which will add additional scatter and uncertainty to the 
relation presented in Fig. 12. 
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FIG. 14. Similar to Fig. 12, but now showing the rescaled 
skewness with /3 — 0.7 plotted against the sum of the neutrino 
masses Sm„. We assume A|j = 2.46 x 1(T 9 (its WMAP5 
value). There is significantly reduced scatter between the 
pressure profiles for this statistic, and it retains a leading- 
order quadratic dependence on Em„. 
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FIG. 15. The rescaled kurtosis {r 4 )/(T 2 > 19 plotted against 
as- It is evident that this statistic is nearly independent of as, 
as expected based on the scalings in Tables I and II. However, 
it is still dependent on the ICM astrophysics, as represented 
by the different pressure profiles. Thus, this statistic — like 
the rescaled skewness — can be used for determining the cor- 
rect gas physics model. 



Nonetheless, we are encouraged by the mild sensitivity 
of S , 3 i ^=o.7 to the complicated astrophysics of the ICM. 
Moreover, it retains a strong dependence on erg. In par- 
ticular, using the results presented in Table I for (T 2 ) 
and (T 3 ), we find that £3,^=0.7 oc Cg~ 6 . For a hxed value 
of erg = 0.817, the scatter in §3^—0.7 between the dif- 
ferent profiles is only 7%. This scatter is much smaller 
than that in the variance or skewness: for the same value 
of cr 8 = 0.817, the scatter in (T 2 ) is ps 50% between 
these profiles, while for (T 3 ) it is « 35%. This statis- 
tic could thus be used for a high-precision determination 
of erg from tSZ measurements, with significantly reduced 
sensitivity to theoretical systematic uncertainties due to 
unknown physics in the ICM. For a full-sky, CV-limited 
experiment, we estimate a SNR « 40 for S3,p=o.7, which 
implies a < 1% constraint on erg, in agreement with the 
estimates presented for £3,0=1.4 earlier. 

However, before applying this technique to real data, 
its insensitivity to the gas physics prescription must be 
tested with detailed numerical tSZ simulations, as we lack 
a good understanding of the theoretical origin of £3,^=0.7- 
The rescaled skewness with j3 = 1.4, though, is well- 
motivated theoretically and provides a clear method to 
constrain the ICM astrophysics in low-mass halos, which 
will greatly reduce the theoretical systematic uncertainty 
in tSZ measurements. We defer an application of the 
P = 0.7 approach to future work. 



Appendix B: Higher-Order tSZ Moments 

In the main text of the paper, we focus on the second 
(variance) and third (unnormalized skewness) moments 
of the tSZ signal. In this appendix, we provide additional 
calculations of the fourth moment (the unnormalized kur- 
tosis) and investigate a "rescaled kurtosis" statistic sim- 
ilar to the rescaled skewness defined in Eq. (|6|). The 
unnormalized kurtosis is simply given by Eq. ^ with 
N = 4. As for the variance and skewness, we compute 
this quantity for each of the various pressure profiles de- 
scribed in the text, and investigate the dependence on 
erg. We define an amplitude A4 and scaling 0:4 precisely 
analogous to A 2 ,3 and 02,3 in Eq. (|5|). The results for 
each of the profiles are given in Table II. 

As for the variance and skewness, the profiles that 
account for feedback processes (Battaglia and Arnaud) 
have a steeper scaling with eg than those that do not, 
because the kurtosis signal for these profiles is dominated 
to an even greater extent by the most massive, rare halos. 
Additionally, note that the scatter between the various 
profiles is even smaller than that for the skewness (and 
much smaller than that for the variance). 

Furthermore, it is possible to construct a "rescaled kur- 
tosis" statistic similar to the rescaled skewness defined in 
Eq. ([5]) , which should be nearly independent of the back- 
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Profile 


Ai [10 s a*K 4 ] 


Qf4 


Arnaud 


4.22 


14.7 


Battaglia 


3.58 


14.5 


Batt. Adiabatic 


5.48 


12.5 


Komatsu-Seljak 


4.70 


13.5 



TABLE II. Amplitudes and power-law scalings with as for the 
tSZ kurtosis. The first column lists the pressure profile used in 
the calculation (note that all calculations use the Tinker mass 
function). The amplitudes are specified at as = 0.817, the 
WMAP5 maximum-likelihood value. All results are computed 
sbv= 150 GHz. 
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FIG. 16. Fraction of the rescaled kurtosis {T 4 >/(T 2 > 1S5 con- 
tributed by clusters with virial mass M < M max . Although a 
non-trivial fraction of the signal comes from low-mass objects, 
comparison with Fig. 6 indicates that the rescaled kurtosis is 
mostly sourced by more massive halos than those responsible 
for the rescaled skewness. 



ground cosmology. In particular, based on the scalings in 
Tables I and II, one can immediately see that the quan- 
tity (T 4 )/(r 2 ) 1 ' 9 effectively cancels the dependence on 
as for each of the profiles. This result is clearly seen 
in Fig. 15, which presents this statistic as a function of 
as- Encouragingly, a clear dependence on the ICM as- 
trophysics persists, as for the rescaled skewness in Fig. 5. 
Note that the same hierarchy is seen between the differ- 
ent profiles in both figures — again, this arises because of 



the suppression of the tSZ signal in low-mass halos for the 
profiles that include significant feedback, which decreases 
their value of the variance, thus increasing their value of 
the rescaled skewness and kurtosis. In addition, it ap- 
pears that the rescaled kurtosis (Fig. 15) is somewhat 
less dependent on erg than the rescaled skewness (Fig. 5), 
although the difference is not very significant over the 
feasible range of as ■ Physically, this is likely due to the 
fact that the kurtosis is dominated even more than the 
skewness by very massive halos, for which the different 
pressure profiles do not differ widely in their predictions. 
Thus, the various profiles have a similar relative scal- 
ing between the kurtosis and the variance, making the 
rescaled kurtosis nearly independent of as for the same 
choice of "rcscaling exponent" . We verify this interpre- 
tation in Fig. 16, which shows the fraction of this sig- 
nal contributed by clusters with virial mass M < Af max . 
Compared to the rescaled skewness (Fig. 6), the rescaled 
kurtosis is sourced more predominantly by fairly massive 
halos. The rescaled kurtosis thus also provides a route 
to determining the gas physics of the ICM, nearly inde- 
pendent of the background cosmology. 

Unfortunately, a measurement of the tSZ kurtosis ap- 
pears extremely challenging with current data, both be- 
cause of confusion with other signals (e.g., the kinetic SZ 
effect, point sources, and so on) and because the cosmic 
variance for this statistic is quite large, as it is sourced 
by very massive clusters. Previous work on the kinetic 
SZ kurtosis allows some rough estimates of this contam- 
ination to be made, however [El|. In [5(j, estimates 
are given of the normalized (dimcnsionlcss) skewness and 
kurtosis for the tSZ and kSZ effects, as determined from 
hydrodynamical simulations. Based on their results, it 
appears that the dimensionless kurtosis due to the tSZ 
effect is at least a few times larger than that due to the 
kSZ effect. Thus, assuming that the variance of the kSZ 
signal is no larger than that of the tSZ signal, this implies 
that the value of (T 4 ) that we have calculated for the tSZ 
effect will be at least a few times larger than the contri- 
bution due to the kSZ effect (at 150 GHz). Hence, the 
kSZ contamination may not be a major problem for a de- 
tection of the tSZ kurtosis. However, contamination from 
point sources will likely remain a significant problem for 
an experiment with only one or two frequencies. Thus, 
ACT and SPT are unlikely to make a detection, though 
Planck may do so. Neglecting any potential contami- 
nants and using the Schgal simulation analysis described 
earlier, we find that a cosmic variance-limited full-sky ex- 
periment could achieve a SNR « 18 for the tSZ kurtosis 
or a SNR « 22 for the rescaled kurtosis. The results 
from Planck will clearly have a somewhat lower SNR 
than these estimates, but it may still be worthwhile to 
investigate these quantities using the forthcoming Planck 
sky maps. 



[1] Sunyaev, R. A., & Zeldovich, Y. B. 1970, Ap&SS, 7, 3 



[2] Reese, E. D., Mroczkowski, T., Menanteau, F., et al. 



14 



2011, arXiv:1108.3343 
[3] Plagge, T. J., Marrone, D. P., Abdulla, Z., et al. 2012, [28 
arXiv:1203.2175 

[4] Lancaster, K., Birkinshaw, M., Gawroriski, M. P., et al. [29 

2011, MNRAS, 418, 1441 
[5] AMI Consortium, Shimwell, T. W., Rodriguez- [30 

Gonzalvez, C, et al. 2011, arXiv:1101.5590 
[6] Marriage, T. A., Acquaviva, V., Ade, P. A. R., et al. [31 

2011, ApJ, 737, 61 
[7] Williamson, R., Benson, B. A., High, F. W., et al. 2011, [32 
ApJ, 738, 139 [33 
[8] Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 

2011, A&A, 536, A8 [34 
[9] Dunkley, J., Hlozek, R., Sievers, J., et al. 2011, ApJ, 739, 

52 [35 
[10] Reichardt, C. L., Shaw, L., Zahn, O., et al. 2011, 

arXiv:1111.0932 [36 
[11] Komatsu, E., & Seljak, U. 2002, MNRAS, 336, 1256 [37 
[12] Cooray, A., Hu, W., & Tegmark, M. 2000, ApJ, 540, 1 
[13] Rubifio-Martm, J. A., & Sunyaev, R. A. 2003, MNRAS, [38 
344, 1155 

[14] Holder, G. P., McCarthy, I. G., & Babul, A. 2007, MN- [39 

RAS, 382, 1697 [40 
[15] Munshi, D., Joudaki, S., Smidt, J., & Coles, P. 2011, 

arXiv:1106.0766 [41 
[16] Fowler, J. W., Niemack, M. D., Dicker, S. R., et al. 2007, 

Applied Optics, 46, 3444 [42 
[17] Swetz, D. S., Ade, P. A. R., Amiri, M., et al. 2011, ApJS, 

194, 41 [43 
[18] Carlstrom, J. E., Ade, P. A. R., Aird, K. A., et al. 2011, 

PASP, 123, 568 [44 
[19] Schaffer, K. K., Crawford, T. M., Aird, K. A., et al. 2011, 

ApJ, 743, 90 [45 
[20] Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, 

ApJS, 180, 330 [46 
[21] Nozawa, S., Itoh, N., Suda, Y., & Ohhata, Y. 2006, 

Nuovo Cimento B Serie, 121, 487 [47 
[22] Cole, S., & Kaiser, N. 1988, MNRAS, 233, 637 
[23] Komatsu, E., & Kitayama, T. 1999, ApJ, 526, LI [ 
[24] Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 

2010, ApJ, 724, 878 [49 
[25] Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80 
[26] Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, [50 

ApJ, 490, 493 [51 
[27] Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vecchia, C. 



2008, MNRAS, 390, L64 

Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 
688, 709 

Battaglia, N., Bond, J. R., Pfrommer, C, & Sievers, J. L. 
2011, arXiv:1109.3711 

Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A 
& A, 517, A92 

Battaglia, N., Bond, J. R., Pfrommer, C, Sievers, J. L., 
& Sijacki, D. 2010, ApJ, 725, 91 
Piffaretti, R., & Valdarnini, R. 2008, A&A, 491, 71 
Shaw, L. D., Nagai, D., Bhattacharya, S., & Lau, E. T. 

2010, ApJ, 725, 1452 

Bhattacharya, S., Nagai, D., Shaw, L., Crawford, T., & 

Holder, G. P. 2012, arXiv: 1203.6368 

Sun, M., Sehgal, N., Voit, G. M., et al. 2011, ApJ, 727, 

L49 

Sehgal, N., Bode, P., Das, S., et al. 2010, ApJ, 709, 920 
Battaglia, N., Bond, J. R., Pfrommer, C, & Sievers, J. L. 

2011, arXiv:1109.3709 

Wilson, M. J., Sherwin, B. D., Hill, J. C, et al. 2012, 
arXiv: 1203.6633 

Ichiki, K., & Takada, M. 2011, arXiv:1108.4688 
Bhattacharya, S., Heitmann, K., White, M., et al. 2011, 
ApJ, 732, 122 

Brandbyge, J., Hannestad, S., Haugb0lle, T., & Wong, 
Y. Y. Y. 2010, J. Cosmology Astropart. Phys., 9, 14 
Marulli, F., Carbone, C, Viel, M., Moscardini, L., & 
Cimatti, A. 2011, MNRAS, 418, 346 
LoVerde, M., & Smith, K. M. 2011, J. Cosmology As- 
tropart. Phys., 8, 3 

Wagner, C, Verde, L., & Boubekeur, L. 2010, J. Cosmol- 
ogy Astropart. Phys., 10, 22 

Dalai, N., Dore, O., Huterer, D., & Shirokov, A. 2008, 
Phys. Rev. D, 77, 123514 

Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, 
ApJS, 192, 18 

LoVerde, M., Miller, A., Shandera, S., & Verde, L. 2008, 

J. Cosmology Astropart. Phys., 4, 14 

Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 

2009, ApJ, 692, 1060 

Planck Blue Book. http://www.rssd.esa.int/SA/PLANCK/ 
docs/Bluebook-ESA-SCI%282005%291.V2.pdf 
Zhang, P., Pen, U.-L., & Wang, B. 2002, ApJ, 577, 555 
Castro, P. G. 2004, Phys. Rev. D, 70, 049902 



